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Abstract 

Using He@C§o as an example, we demonstrate that static potential of the fullerene core es- 
sentially alters the cross section of the two-electron ionization differential in one-electron energy 
dcr ++ {uj) / de . We found that at high photon energy prominent oscillations appear in it due to 
reflection of the second, slow electron wave on the Cqq shell, which "dies out" at relatively high 
e values, of about 2-^3 two-electron ionization potentials. The results were presented for ratios 
Rc ao (uj,e) = da ++ (u, e)/da a++ (oj, e), where da a++ {uj , e) / de is the two-electron differential pho- 
toionization cross section. We have calculated the ratio Rij u i = af + {uj)/a^ ++ {uj), that accounts 
for reflection of both photoelectrons by the Cqq shell. We have calculated also the value of two- 
electron photoionization cross section <j ++ (uj) and found that this value is close to that of an 
isolated He atom. 

PACS numbers: 32.80.Fb, 31.15.Ja 
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I. INTRODUCTION 



Elimination of two electrons by a single photon from an atom or multi-atomic formation, 
such as cluster or fullerene, can take place only if the inter-electron interaction is taken into 
account. 

The desire to study the manifestation of this interaction stimulates extensive experimental 
and theoretical investigation of the process that has particularly intensified during last ten- 
fifteen years (see, e.g., Although a number of atoms were studied, the primary 
attention was given to Helium. At this moment theoretical and experimental investigations 
cover the frequency region from the near threshold region up to photon energy uj much 
higher than the two-electron ionization potential I ++ . 

Very close to threshold the two-electron photoionization cross section is determined by 
so-called Wannier regime jjj with both electrons strongly repulsing each other that results 
in acquiring almost the same energy and moving in opposite directions. 

With increase of u for uj considerably more than I ++ , so-called shake-off (SO) mechanism 
became dominant, in which one electron leaves the atom carrying away almost all energy 
uj, while the second is removed due to alteration of the field acting upon it after the first 

n 

electron emission. With energy growth so called quasi-free (QF) mechanism |5( becomes 
increasingly important. This mechanism accounts for almost equal sharing of photon en- 
ergy between photoelectrons, where their interaction on the way out of the ionized atom is 
inessential. 

Obviously, the cross section of two-electron photoionization <j ++ (uj) is determined by 
initial and final state wave function of the considered object - target atom at the beginning 
and ion with two continuous spectrum electrons at the end. The situation is simplified 
considerably at high enough uj, where cross sections and other characteristics of the two- 
electron ionization are expressed via the initial state wave function only. 

In fact, one has to have in mind that two-electron ionization is a pure three-body problem 
only for two-electron atom and ions. This process for any other, more complicated objects 
is determined by more sophisticated wave functions. This is why most of the attention is 
given usually to studies of two-electron photoionization of H~ , He and helium-like atoms. 
At high uj the a ++ and a + have similar u dependence j^J. Therefore it is convenient to 
characterize the process by the ratio R(uj) = a ++ (uj) / a + (uj) . 
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The SO leads for uj » I ++ to R(u) = Rso (see |?J and references therein). Inclusion 
of QF increases R(uj) considerably but at high enough uj the ratio again reaches its uj- 
independent value, Rqf > Rso 



During last time a lot of attention is given to photoionization o 



not isolated atoms, but 



!. The research in 



atoms encapsulated into the fullerene, mainly C 60 shell (see, e.g. 
this area until now is pure theoretical. But we are positive, that it will become an object of 
experimental studies in not too distant future. It is well known, that Cqq radius R is much 
bigger than that of the atoms staffed inside the so-called endohedral atom, A@C 60 . 

The inclusion of the Cqo shell can affect the photoionization atom in four directions. At 
first, and this is usually taken into account, the photoelectron emitted by A is reflected by 
the C*6o shell that "works" as a static potential. At second, the Cqq shell has its internal 
degrees of freedom [11] that can be excited thus modifying the photoionization of the atomic 
shells in A@Cqq in a similar way that one shell affects the photoionization of the other in 
isolated atoms jlj ]. Among the Cqo degrees of freedom the most prominent is the giant 
resonance that persists not only in Cqq jllj], but in its ions as well Q|. Real or virtual 
excitations of the giant resonance can modify the A atom photoionization cross section from 
almost complete screening at uj — » to prominent enhancement, at uj in the vicinity of 
the giant resonance maxima 14, Q|. Note, that at high uj this screening becomes mainly 



inessential. 

Other two directions of Cqq influence upon the A atom photoionization is the direct knock 
out of Cgo electrons by the photoelectron from A and the participation of Cqo electrons in 
the decay, both radiative and non-radiative, of the vacancy, created in A after photon 
absorption. Note that this mechanism is similar to one of the frequently discussed channels 
of the molecular vacancy decay. The corresponding two possibilities in application to Cqq 
are not yet studied at all. 

In this paper we consider a more complex process, namely two-electron photoionization 
of He@CQQ. We will concentrate on the high uj region. Of course, to investigate this process 
experimentally is more complicated than one-electron photoionization. However, there are 
no indications whatsoever that such an investigation is either impossible (or even extra- 
ordinary difficult) or uninteresting. 

The role of Cqq shell in two-electron photoionization is more complex than in the single- 
electron, since outgoing electrons can be either both from the inner atom A, or both from the 
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6*60, or one from A and the other from Cqq- However, by measuring two outgoing electrons 
in coincidence we can distinguish all the different processes. 

We will consider mainly high uj. Having in mind that our aim is study the C 60 influence 
upon double ionization, we have to consider the SO mechanism only. Indeed, the QF leads 
to two fast electrons, which are not affected by Cqo, while in SO one of the electrons is slow. 
Its probability to leave He@Ceo can be therefore strongly affected. 

Our attention in this paper will be given to differential in energy cross sections 
da ++ (uj)/de s , where e s is the "slow" photoelectrons energy. It is implied that the "fast" 
electron energy £/ is given by the conservation law ef = uj — e s — I ++ . Atomic system of 
units is accepted in this paper: e = m e — H — 1, where e and m are electron charge and 
mass, respectively. 

We will assume that both electrons are removed from the He atom, thus intensionally 
neglecting a process that one can call " A@C 60 shake off', in which after "fast" electron leaves 
A, the "slow" electron is emitted by the C 60 shell instead of located inside atom A. The 
" A@C eo shake off' is potentially very important due to big number of available for removing 
electrons. Therefore, it is possible that " He@Ceo shake off' is much more probable than the 
ordinary one. It is possible, that an important role can be played by multiple " A@Cqq shake 
of, in which not one but several electrons can be removed from Ceo after instant creation 
of a vacancy in A. One has to have in mind that for an isolated atom A the photoionization 
with excitation cross section <j + *(uj) can be of the order or even bigger than <j ++ (uo). It is 
essential to have in mind that the cr + *(uo) for A@C eo can be converted into a ++ (uj), since 
the excitation energy of A + , particularly in the case of He, is bigger than the ionization 
potential of C 60 . We plan, however, to concentrate on " A@C 60 shake off' in another paper. 

II. MAIN FORMULAS 

The two-electron photoionization cross section of an atom A in initial state i, a1 ++ (uj), 
can be presented by the following expressions: 



where at high uj, uj » I ++ , the differential in energy cross section da^ ++ (uj, e)/de is 




(1) 
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given by the expression 16]: 



daf ++ {uj,e) 32^2 Z 2 ti . 



J *i(0,r)^ (r)dr| 2 . (2) 



de 3coo 7 / 2 

Here ^i(0, r) is the ri = value of ^(ri, r2), which is in our case the initial state wave 
function of atomic He. The " slow" outgoing electron is described by a pure Coulomb wave 
function that describes an s -wave electron that moves in the field of a nucleus with charge 
Z. 

Similar to Eq.Q the cross-section of ionization with excitation <j + *(oj) is given by the 
following expression: 



daf ++ (oo,e) _ 32^2Z 



2 r 2 

2^1 / #i(0,rKo(r)dr| . (3) 



de 3cuj 7 / 2 

n>0 

Here n is the principal quantum number of the second electron excitation level. 

Now let us consider He@Cgo two-electron photoionization and take into account the Cqq 
shell. It is clear that Cqq potential does not affect the "fast" electron, which remains to be 
described by a plane wave. The C§o shell generates potential field that is almost uniform at 
atomic distances r a , where the initial state of He is located. Therefore, for the initial state 
the embedding of an atom A into Cqq can lead only to the shift of the energy scale, i.e., to 
modification of the ionization threshold. This feature is inessential for us, since uj » I ++ . 

Thus, we have to take into account the action of C^o upon the slow outgoing electron. The 
fullerene is a very complex structure. Therefore, to take into account its action upon photo- 
electron essential simplifications are necessary. We will follow here the approach developed 
in a number of papers (see 0,0] and references therein), that substitutes the complicated 
C 60 multi-atomic structure by a very simple so-called bubble potential V(r) = V 6(r — R), 
where Vq is chosen to reproduce the affinity energy of the negative ion Cqq, and 5(r — R) is 
Dirac delta function. 

As it was demonstrated in Q and Q in connection to the one-electron photoionization, 
the influence of the bubble potential can be taken into account analytically, by constructing 
the outgoing wave function as a superposition of regular y?d(r) and irregular Xei( r ) (singular 
at r —>■ 0) solutions of the Schrodinger equation for an electron with energy e. Inside 
the potential bubble the photoelectron wave function ip e i(r) differs from <f e i(r) only by a 
normalization factor Di(e) : ip £ i{ r ) — Di(e)(p £ i(r). The factor Di(e) depends upon the 
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photoelectron energy e, with I being its electron angular momentum. Outside the ^-sphere 
the function ip £ i(r) is a linear combination of y?d( r ) and Xd( r )- The coefficients of the linear 
combination are defined by the matching conditions of the wave functions on the spherical 
shell, i.e., at r = R. 

In the (^-potential approximation the differential two-electron photoionization of He@Cso 
is thus given by formula similar to that includes, however the factor Di(e): 



7T 



D (e)\ 2 \ j ^(0,r)^ (r)rfr| 2 . (4) 



de 3cu) 7 / 2 
The expression for |.D;(e)| 2 is derived in |l9l. l2flj] and is presented as 

| DM | 2 = (fe/AL) 2 () 

Here u k i(R) and Vki(R) are functions connected to the radial parts of regular (p £ i(r) and 
irregular Xei( T ) functions by relations Uki(r) = r^ ei (r) and ^jy(r) = rx £ i{r); k = v2e, AL is 
the discontinuity of the logarithmic derivative of the wave function at r = R, connected to 
the fullerene radius R and the electron affinity If of the empty Cqq through the expression 

AL = -(3(1 + coth/3ii!), (6) 

where [3 = \/2If- The formula obtained are valid for low enough energies of the "slow" 
photoelectron. Namely, its wave length should be much bigger than the thickness of the 
fullerene shell. 

At first one should calculate the ratio of differential cross sections (jJJ) to (J2J) Rc 60 ( £ ) that 
according to Eqs.® and (@J is independent upon u 

R C60 (e) = dat + (u,e)/da« ++ (u,e) = \D (e)\\ (7) 

and then turn to the ratio Rij u i of cross sections <j ++ (uo) and a a++ (uo) that are determined 
using Eq.(Ul with daf + (u,e)/ds given by Eq.(JH) and Eq.(J2J), respectively. 

R lJul = at + (co)/at ++ (co) (8) 
It is evident from (J7J) that determined in this way Rij u i is u -independent at high ui. 
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Since <p e o( r ) and Xei( r ) are pure Coulomb functions, the functions Uki and Vki (from 
Eqs. (|6l7jl ) are expressed via the regular Fi(rj,p) and irregular Gi(rj,p) Coulomb wave func- 
tion, respectively. Namely, one has w^(r) = F^—j, kr) and VbAr) = Gi(—j, kr), i.e. regular 
and irregular Coulomb wave functions that can be found in [20]. 

Let us note that the formula (jZJ can be extended to lower u. Indeed, the interaction 
between two ionized electrons is particularly essential near ionization threshold to > I ++ , 
where Wannier regime prevails. It is known that already several eV above threshold the 
Wannier expressions are no more valid. Therefore, it is reasonable to assume that at about 
lOeV the interaction between outgoing electrons is inessential. One of them is represented 
by an s-wave, while the other by a p-wave, in order to preserve the total initial state angular 
momentum that is zero for He@CQo and I = 1 for the photon. Due to action of the Ceo shell 
the second electrons wave function is modified in the same way as y? £ o( r ) i n (2). Thus, it is 
reasonable to expect that to take this action into account one has to introduce the factor 
| Z^i (^) | 2 into ((Zj). As a result, one obtains instead of Eq.fllJ) the following expression: 

Rcfae) = da+ + (uj } e)/da? ++ (uj 1 e) = |£> (e)| 2 |I>i(a; - / ++ - e)\ 2 . (9) 

Since for low and medium energies the cross section da1 ++ (u,e)/de cannot be calculated 
using Eq.(J2J), we have to use other sources of absolute values of it in order to obtain Rc 60 (uj). 

As such, either experimental data or results of existing calculations of differential in 
energy e cross section can be used to substitute into 

R OM (cj)=r I |Do(^)| 2 |A(^-/ ++ -^)| 2 rf< ++ (^^)/ f 1 dat ++ (u,e). (10) 

Jo Jo 

In order to have a crude estimation of the role of both factors in @, we have calculated 
in Eq. (fTU|) the factor Rc 60 (uj) using Eq.((TJ) well outside its range of validity uj 3> which 
gives us R%1(uj). 

Since at u close to daf ++ (cu, e)/de is almost cu-independent, as an estimation for 
Rc 60 (uj) at uj > I ++ can serve an approximate relation 

4LH = 1 \Da(£)\ 2 \Dx(" ~ I ++ ~ e)\ 2 de/{uj - (11) 

J o 
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III. DETAILS OF CALCULATIONS 



First of all, one should notice that the term k/AL in denominator of Eq.© was taken 

nn 

with the sign "minus", whereas there was "plus" in the Ref . |18L [19|. This is because the 
authors of the latter references used the irregular Coulomb functions vui(r) = —Vki(r), whilst 
we used the Coulomb functions Ujy(r) presented in the Handbook [20J. 

The Eqs.fCTTTj) include the function Di(e) at e = 0, for e = uj — I ++ . It is clear that 
the Coulomb functions Fi(—j, kr) and Gi(—j, kr), presented in the expression (0) for Df, 
have singularities at zero energy, i.e. at k = that corresponds to e = 0. However, 
the presentations of Coulomb functions given in [2£j] enable to obtain the following limit 
relations: 



lim^ = fi(Z,r) = s/^J 2l+1 (V8Z^), (12) 
fc-o \/k 

lira ^ = 9o(Z, r) = - Jl [2 7 - 1 + In(2Zr)] Jx(Vs^). (13) 

fc ^o Vk V2tcZ V 7T 

Here 7 is Euler's constant, J M (x) is a usual Bessel function of the first kind, and p(a) is the 
following sum, which is infinite but rapidly convergent: p(a) = Yl°iLoPj{ a )i where Pj(ct) 
could be obtained by the recurrence relations 



P 



, 1) (j!)2 +^-x 



(14) 



J'(J " 1) 

with the initial values Po = 1, -Pi = 0. Functions gi(Z,r) for I > could be then obtained 
by using Eqs. (jl2j) - (jl4j) and the following recurrence relation 

g l {Z,r)f l ^ 1 {Z,r) = g^Zrfffar) + l -. (15) 



The latter is a consequence of the Wronskian relation (see Eq. (14.2.5) in Ref. 20]) for the 
regular and irregular Coulomb functions. Note that Eq. (jl2j) for the regular Coulomb func- 
tions could be found in numerous handbooks on non-relativistic quantum mechanics (see, 
e.g. Ref.([13|)). within the accuracy of a normalization factor. However, we did not find in 
scientific literature the corresponding relations (|13 )) -()15 )) for irregular Coulomb functions. 

Using the limit relations (|12j) - ([TK|) , one could present Eq.(JHJ for the zero energy in the 
form: 

A2(0) = [(AL)fl(Z, R)fi(Z, R) - l] 2 + (AL)*tf(Z, R) ' (16) 



In particular, for the helium atom expression (j!6|) yields: -Dq(O) = 1.38590 and Df(0) = 
3.09129. 

The function can be obtained using variational principle that minimizes the initial 
state energy. One has to have in mind, that, generally speaking, this procedure is able to 
reproduce initial state energy absolutely accurate, leading to wave function that satisfies the 
precise Schrodinger equation for two electrons in the field of a nucleus with the charge Z. In 
its practical implementations, however, the variation wave functions reproduces the initial 
energy approximately and therefore is good on the average. 

As to photoionization, since one of the outgoing electrons is fast, its wave function is a 
plane wave with a wavelength short in the atomic scale. Therefore, photoeffect is able to 
test the short range behavior of ^(r^, r 2 ), namely, as it is seen from Eq.(Q), ^AO, r). 

tion I2J, 



23, 



24, 



This is why instead of a variational, we use a locally correct wave functi 
that describes with high accuracy the two-electron photoionization [16( and some average 
characteristics of the initial state. 



IV. RESULTS OF CALCULATIONS 

Here we presented the results of our calculations for the two-electron photoionization of 
He@CQo- For this object the radius R = 6.64a.u. and the electron affinity of the empty 
C*6o is equal to If = 2.65eV. The initial state of He considered in this paper is i = Is, i.e. 
I = 0. The results for Rc 60 (e) from Eq.® are presented in Fig.l together with the function 
D\(e) from the Eqs.(j9HTTj). As it was expected and similar to the case of one-electron 
photoionization, the ratio Rc 60 ( e ) ls a strongly oscillating function of e. 

Fig. (2) depicts Rc 60 (u,e), given by Eq.fjH} for several uj as functions of e. It is not 
incidental, that the curves behavior at e = and e = uj — 1 ++ are different. This is because 
one of the electron is represented by an s-wave while the other by p-wave. 

Note that Rc 60 {uj, e) becomes a prominently variating curve already for uj — 1 ++ > 20eV. 
For u — I ++ > QQeV it has many oscillations. 

Fig. (3) presents the results for the ratio Rij u i{uj) determined by Eq.(JHJ). As one could 
see, with uj growth this ratio rapidly approaches the asymptotic value Rij u i — > 1. Note, 
that in Eq.fjHJ) we assumed that "slow" electron wave function is modified by the Cqo shell. 
In fact, this is incorrect already for uj > lQOeV. 
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FIG. 1: The reflection and refraction factors D 2 L [L = 0, 1) as a function of photoelectron energy 
e for He@C 60 . 

Fig. (4) presents the relations Rc 60 (u>) given by Eqs. (jlOMfTj) for both e distributions, that 
are valid respectively at high, and at low uj. It is seen that both curves, Rq &q (u) and Rq so (u), 
prominently oscillate up to u ~ lSOeV^ and are close to each other. It permits us to suggest 
that the curve Rc 60 (uj) is valid for correct e distribution instead of only their limiting cases 
for high and low uj. The strongest is the role of the fullerene shell at u < HOeV. 
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FIG. 2: The cross sections ratio Rc 60 {^,£) as a function of photoelectron energy e for different 
values of photon energy u. 
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FIG. 3: The cross sections ratio Rijul as a function of photon energy to. 
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FIG. 4: The cross sections ratios Rq^ and Rq as functions of oj for high and low photon energy, 
respectively. 
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